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ABSTRACT 


As  an  example  of  ships  used  as  weapons  (SAW),  an  oil  tanker  is  hijacked  and 
commandeered  by  terrorists  to  collide  with  a  high-value  maritime  or  shore  target.  If  sunk 
or  destroyed  in  a  shipping  lane  as  a  result  of  a  counter  measure,  the  SAW’s  collateral 
damage  would  severely  disrupt  the  traffic  flow  in  the  shipping  lane.  To  prevent  such  a 
disruptive  catastrophe,  non-destructive  measures  must  be  implemented  to  cause  the  SAW 
to  deviate  from  its  destructive  path  toward  the  target.  One  such  a  measure  involves  a 
strategic  application  of  forces  induced  by  water  plume  barriers  (WPB)  to  the  SAW.  The 
goal  of  this  thesis  is  to  examine  the  feasibility  of  realizing  such  a  measure. 

Toward  this  goal,  a  mission  analysis,  using  the  Singapore  Strait  as  setting  and 
petrochemical  plants  on  Jurong  Island  as  targets  of  a  SAW  attack,  establishes  the 
requirement  on  the  deviation  of  the  SAW  path  from  its  destructive  course.  The  nominal 
WPB-induced  force  that  satisfies  the  deviation  requirement  is  estimated  using  ship 
hydrostatics.  Solving  the  equations  of  motion  governing  the  response  of  the  SAW  to  a 
strategic  application  of  a  WPB-induced  force  yields  the  SAW’s  motion,  which  is  used  to 
define  a  range  of  the  WPB-induced  forces  and  their  application  locations  and  durations 
that  satisfy  the  SAW’s  path  deviation  requirement. 

Parametric  studies  were  conducted  for  a  range  of  physically  realizable  WPB- 
induced  forces  and  application  times.  The  results  demonstrate  that,  in  principle,  the 
objectives  of  this  work  are  achievable.  These  results  will  be  validated  upon  the 
completion  of  an  on-going  research  by  National  University  of  Singapore.  The  range  of 
the  WPB-generated  forces  and  the  application  durations  serve  as  requirements  to  the 
generation  of  water  plume  barriers. 
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I.  INTRODUCTION 


A.  BACKGROUND 

This  chapter  discusses  the  increasing  concerns  over  maritime  terrorist  attacks,  the 
research  motivation,  and  the  objective  of  this  thesis.  It  also  defines  the  research  scope, 
formulates  the  problem  to  solve  in  this  research,  discusses  the  approach  to  solving  the 
problem,  and  outlines  the  structure  of  the  thesis. 

“A  major  terrorist  attack  that  closed  a  port  ...  for  weeks  would  have 
severe  economic  consequences  on  world  trade  because  it  would  inflict 
major  disruptions  in  complex  just-in-time  supply  chains  that  comprise  the 
global  economy,”  the  World  Economic  Forum  said  in  its  Global  Risks 
2010  report,  released  in  January  [1]. 

U.S.  Congressional  Research  Service  has  reported  to  U.S.  Congress  that  the  threat 
of  maritime  terrorism  is  significant  and  takes  many  forms  of  variation  and  raised 
concerns  about  the  possibility  of  maritime  terrorist  attacks  [2],  Maritime  vessels  and 
facilities  are  particularly  vulnerable  to  terrorism.  Several  planned  seaborne  attacks 
occurred  in  the  past  decade.  Two  notable  events  are  the  explosives-laden  dinghy  attack 
against  the  USS  Cole  in  2000,  resulting  in  17  casualties,  and  an  attack  on  the  Limburg  in 
2002,  killing  one  crewman  and  spilling  90,000  barrels  of  oil  into  the  Gulf  of  Aden  [1]. 

Maritime  terrorist  attacks  are  not  just  confined  to  any  specific  region.  The  Straits 
of  Malacca  and  Singapore  between  peninsular  Malaysia  and  Sumatra  can  potentially 
attract  such  attacks  [3],  One  of  the  world’s  busiest  shipping  lanes,  the  Straits  were  used 
by  more  than  70,000  ships  in  2007.  Up  to  80%  of  China’s  oil  imports  and  90%  of  Japan’s 
crude  oil  imports  pass  through  the  Straits  [1].  They  are  thus  vital  passageways  for  the 
transportation  of  petroleum  products  from  the  Middle  East  and  East  Asia.  The  U.S. 
Energy  Infonnation  Administration  has  identified  the  Straits  as  one  of  the  world’s  most 
strategic  world  oil  transit  chokepoints  (critical  parts  of  global  energy  security)  [4]. 

If  terrorists  hijacked  an  oil  tanker  and  used  it  as  a  weapon,  hence  ship  as  a  weapon 
(SAW),  in  a  successful  maritime  attack  in  the  Straits,  its  collateral  damage  would  disrupt 
the  shipping  lanes  and  the  petrochemical  sector  in  particular,  thereby  bringing  devastation 
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to  not  only  the  local  economies  but  also  the  global  economy  as  well.  A  defensive  measure 
is  thus  needed  to  stop  the  terrorists  from  achieving  their  mission.  Such  a  defensive 
measure  to  counter  the  SAW,  however,  must  not  sink  or  destroy  the  SAW  in  a  shipping 
lane,  because  its  collateral  damage  would  severely  disrupt  the  traffic  flow  in  the  shipping 
lane.  To  prevent  such  a  disruptive  catastrophe,  non-destructive  measures  must  be 
implemented  to  cause  the  SAW  to  deviate  from  its  destructive  path  toward  the  target. 
One  such  a  measure  involves  a  strategic  application  of  forces  induced  by  water  plume 
barrier  (WPB)  to  the  SAW. 

A  WPB  is  the  free  water  surface  that  is  being  pushed  up  by  the  bubbles  created 
by,  for  example,  underwater  explosions.  These  bubbles  are  strategically  located  to  create 
a  WPB  of  a  particular  shape.  The  resulting  WPB  then  exerts  a  force  to  the  SAW  to  cause 
the  SAW  to  move  in  a  direction  away  from  its  intended  destructive  course.  Chapter  III 
discusses  the  formation  of  a  WPB  with  bubbles  and  its  interaction  with  the  hull  of  the 
SAW. 

B.  RESEARCH  MOTIVATION 

The  motivation  of  this  research  is  to  explore  the  feasibility  of  using  a  WPB  as  a 
non-destructive  measure  to  destroy  the  positional  stability  (ship’s  ability  to  return  to  its 
original  path)  of  a  SAW  aimed  towards  a  high-value  target.  Such  a  non-destructive 
measure  should  result  in  neither  sinking  nor  exploding  the  SAW. 

C.  RESEARCH  OBJECTIVE 

The  research  reported  in  this  thesis  is  part  of  a  multi-year  research  project.  This 
research  project  is  funded  by  National  University  of  Singapore  (NUS)  for  a  fund  of 
approximately  SGD  400,000  over  three  years  (2010-2013)  and  the  Naval  Postgraduate 
School  (NPS)  for  a  fund  of  USD  100,000  over  a  year  (2010-2011).  The  ultimate 
objective  of  the  research  project  is  to  develop  and  deploy  a  Water  Plume  Barrier  (WPB) 
system  to  counter  a  SAW. 
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This  thesis  serves  to  initiate  the  research  project;  its  objective  is  to  assess  the  feasibility 
of  applying  WPB-induced  forces  to  the  hull  of  the  SAW  to  alter  its  course  and,  thereby, 
to  generate  the  preliminary  requirements  on  the  magnitude  and  sustainment  of  the  force 
generated  by  a  WPB.  Specifically,  it  attempts  to  solve  the  problem  stated  in  Section  D. 

D.  PROBLEM 

A  SAW  is  commandeered  by  terrorists  to  ram  a  petrochemical  processing  plant 
located  at  the  southernmost  end  of  Jurong  Island  in  the  Singapore  Strait.  If  sunk  or 
destroyed  as  a  result  of  a  counter  measure,  the  SAW’s  collateral  damage  would  severely 
disrupt  the  traffic  flow  in  the  shipping  lane.  To  prevent  such  a  disruptive  catastrophe, 
non-destructive  measures  must  be  implemented  to  cause  the  oil  tanker  to  deviate  from  its 
destructive  path  toward  the  target.  One  such  a  measure  involves  a  strategic  application  of 
WPB-induced  forces  to  the  SAW.  The  feasibility  of  employing  the  WPB  measure  needs 
to  be  assessed.  Assessing  it  amounts  to  answering  the  research  questions  posed  in 
Section  E. 

E.  RESEARCH  QUESTIONS 

The  primary  research  question  is: 

Is  it  feasible  to  use  a  WPB  as  a  means  to  counter  a  SAW? 

Answering  this  question  amounts  to  answering  the  following  questions: 

1)  What  are  the  magnitudes  of  the  WPB-induced  forces  required  to 
destroy  the  SAW’s  positional  stability? 

2)  Where  on  the  hull  of  the  SAW  should  the  WPB-induced  forces  be 
applied? 

3)  For  how  long  does  the  application  of  the  WPB-induced  forces  need  to 
be  maintained? 

F.  RESEARCH  SCOPE 

The  Singapore  Strait  is  the  setting,  and  the  corresponding  operational 

environment  is  considered.  An  oil  tanker  is  used  as  a  SAW  and  the  specific  target  is  a 
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petrochemical  processing  plant  located  at  the  southernmost  end  of  Jurong  Island.  WPB 
optimization  and  simulation  and  WPB  system  design  are  not  covered  in  this  thesis. 


G.  APPROACH  TO  SOLVING  THE  PROBLEM 

A  three-stage  approach  is  employed  to  solve  the  problem  or  to  answer  the 
questions  posed  in  Section  E.  In  the  first  stage,  a  mission  analysis  is  performed.  It 
involves  (1)  an  analysis  of  the  operational  environment  pertaining  to  the  Singapore  Strait, 
(2)  the  postulation  of  a  potential  target,  (3)  the  identification  of  the  threat  (SAW) 
characteristics,  (4)  the  definition  of  a  specific  SAW  scenario,  and  (5)  the  determination  of 
the  SAW’s  altered  path  required  to  achieve  the  counter-SAW  mission.  The  outcomes  of 
(5)  are  criteria  used  in  the  assessment  of  the  feasibility  of  using  a  WPB  to  counter  a 
SAW. 

In  the  second  stage,  a  parametric  study  is  conducted  to  detennine  the  SAW 
motion  resulting  from  responding  to  various  combinations  of  the  three  parameters — the 
magnitude  of  the  WPB-induced  force,  the  location  of  its  application  along  the  hull  of  the 
SAW,  and  the  duration  of  its  application.  This  study  involves  fonnulating  the 
mathematical  problem  of  detennining  the  response  (including  rotation  and  translation)  of 
a  ship  exerted  by  a  WPB-induced  force,  solving  the  resulting  mathematical  problem 
using  MATLAB  to  obtain  the  SAW  motion  and  the  angles  of  deflection  corresponding  to 
various  combinations  of  the  three  parameters.  The  mathematical  formulation  of  the 
response  involves  setting  up  the  equations  of  motion  governing  the  SAW  dynamics, 
defining  initial  conditions,  determining  the  parameters  in  the  equations  of  motion,  such  as 
the  ship  characteristics,  the  magnitude  of  the  applied  WPB-induced  force,  and  the  ship 
hydrodynamics  derivatives.  A  computer  program  is  used  to  generate  the  ship 
hydrodynamics  derivatives,  given  the  SAW  characteristics  and  its  operating  environment. 
Various  magnitudes  of  the  applied  WPB-induced  forces  are  derived  from  a  nominal 
WPB-induced  force  estimated  using  ship  hydrodynamics. 

In  the  final  stage,  based  on  the  parametric  study  results  obtained  in  the  second 
stage,  the  requirements  on  the  magnitude  of  the  WPB-induced  force  and  the  duration  and 
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the  location  of  the  application  of  the  WPB-induced  force  are  established  for  a  successful 
counter-SAW  mission. 

H.  THESIS  STRUCTURE 

This  thesis  consists  of  seven  chapters.  Chapter  I  provides  the  background  of 
maritime  terrorism  as  well  as  the  motivation,  the  objective,  and  the  scope  of  the  research. 
The  problem  and  the  approach  to  solving  the  problem  are  also  described  in  this  chapter. 
Chapter  II  describes  the  mission  analysis  and  its  results,  which  serve  as  the  criteria  used 
in  the  assessment  of  the  feasibility  of  using  a  WPB  to  counter  a  SAW.  Chapter  III 
discusses  the  formation  of  a  WPB  and  its  potential  application  in  a  counter-SAW 
mission. 

Chapter  IV  discusses  the  mathematical  formulation  of  the  ship  response,  which 
involves  setting  up  the  equations  of  motion  governing  the  SAW  dynamics,  defining 
initial  conditions,  detennining  the  parameters  in  the  equations  of  motion,  such  as  the  ship 
characteristics,  the  magnitude  of  the  applied  WPB-induced  force,  and  the  ship 
hydrodynamics  derivatives.  It  also  explains  the  calculation  of  the  nominal  WPB-induced 
force  and  demonstrates  the  solution  implementation  using  MATLAB.  Chapter  V 
discusses  the  parametric  study  and  its  results — the  effects  of  the  magnitudes  of  the  WPB- 
induced  forces  and  the  durations  and  the  locations  of  the  applied  forces  on  the  SAW’s 
orientation  and  kinematics.  The  parametric  study  aims  to  determine  a  combination  or 
combinations  of  these  three  parameters  that  satisfy  the  counter-SAW  mission  success 
requirement. 

Chapter  VI  provides  a  preliminary  assessment  of  the  feasibility  of  using  a  WPB  to 
counter  a  SAW,  in  terms  of  the  magnitude  of  the  WPB-induced  force  and  the  duration 
and  the  location  of  the  force  for  a  counter-SAW  mission  success  requirement.  It 
summarizes  the  results  of  the  different  combinations  of  the  three  parameters  mentioned  in 
Chapter  V. 

Finally,  Chapter  VII  recapitulates  the  problem  and  the  approach  to  solving  the 
problem  and  provides  a  summary  of  the  research  results  and  recommendations  for  future 
research. 
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II.  MISSION  ANALYSIS 


A.  OPERATIONAL  ENVIRONMENT 

The  mission  analysis  in  this  work  consists  of  determining  the  operating 
environment  pertaining  to  the  Singapore  Strait,  considering  a  potential  target,  identifying 
the  threat  (SAW)  characteristics,  defining  a  specific  SAW  scenario,  and  detennining  the 
SAW’s  altered  path  required  to  achieve  the  counter-SAW  mission.  The  SAW’s  altered 
path  serves  a  criterion  used  in  the  assessment  of  the  feasibility  of  using  a  WPB  to  counter 
SAW. 


1.  Geography 

Singapore,  engulfed  by  the  Singapore  Strait,  is  situated  in  the  southernmost  part 
of  peninsular  Malaysia  and  north  of  Indonesian  Riau  Islands  in  South  East  Asia.  As 
Figure  1  shows,  with  the  Malacca  Strait  in  the  west  and  the  South  China  Sea  in  the  east, 
the  Singapore  Strait  is  a  105-kilometer-long,  16-kilometer- wide  strait,  and  it  is  one  of  the 
busiest  waterways  used  by  international  shipping.  The  narrowest  point  of  this  waterway 
is  the  Phillips  Channel,  which  is  only  1 .7  miles  wide  at  its  narrowest  point.  This  creates  a 
natural  bottleneck,  with  the  potential  for  a  collision,  grounding,  or  oil  spill  [4]. 
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Figure  1.  Map  of  Singapore  and  Singapore  Strait  (From  [5]) 

2.  Maritime  Conditions 

Surrounded  by  peninsular  Malaysia  and  the  islands  of  Indonesia,  the  Singapore 
Strait  is  sheltered  from  the  harsh  currents  of  Indian  Ocean  and  South  China  Sea.  The  sea 
condition  in  the  Singapore  Strait  is,  therefore,  relatively  benign  and  the  sea  state  normally 
does  not  exceed  3,  corresponding  to  smooth  wavelets  with  a  wave  height  of  no  more  than 
0.5  meters  [6].  The  mean  surface  wind  speed  can  reach  up  to  2.5  m/s  and  the  average 
non-tidal  streams  are  about  0.5  m/s.  During  monsoon  seasons  the  tidal  streams  can  reach 
up  to  2.5  m/s.  The  surface  water  temperature  is  quite  uniform,  usually  less  than  1°  C  in 
variation  and  the  thermal  effect  of  the  underwater  current  is  thus  minimal  [7]. 

3.  Climate  and  Meteorological  Conditions 

Singapore  is  located  one  degree  north  of  the  Equator.  It  has  a  tropical  rainforest 
climate  with  no  distinctive  seasons,  uniform  temperature  and  pressure,  high  humidity, 
and  abundant  rainfall.  Singapore’s  weather  is  warm  and  humid  all  year  round  and  its 
average  annual  rainfall  is  around  2,340  mm  [8].  However,  it  has  four  periods  of  monsoon 
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seasons — Northeast  Monsoon  Season  (December  to  Early  March),  Inter-monsoon  Period 
(Late  March  to  May),  Southwest  Monsoon  Season  (June  to  September)  and  Inter¬ 
monsoon  Period  (October  to  November)  [9]. 

B.  POTENTIAL  TARGET 

Strategically  located  in  the  Straits  of  Malacca  and  Singapore,  Singapore  is  the 
world’s  busiest  port  and  is  the  largest  transshipment  hub.  The  Singapore  Strait  is  its 
lifeline  for  trade,  food  supply,  and  other  material  needs  [10]. 

Singapore  is  also  one  of  the  world’s  top  petrochemical  hubs.  This  reputation  is 
achieved  with  the  development  of  Jurong  Island  into  a  premier  petrochemical  hub  that 
hosts  over  95  global  companies,  including  Shell,  ExxonMobil,  Chevron,  DuPont™, 
BASF,  Sumitomo  Chemicals,  and  Mitsui  Chemicals.  Jurong  Island  has  drawn 
cumulative  fixed  asset  investments  of  over  S$30  billion  and  employed  about  8,000 
employees  [11].  The  petrochemical  sector  has  contributed  almost  5%  of  Singapore’s 
gross  domestic  product  in  2007  [12].  In  addition,  Singapore  is  the  world’s  top  three  oil 
refining  centers.  Jurong  Island  is  not  just  important  to  Singapore;  it  is  also  critical  to  the 
world  oil  trade. 

The  Straits  of  Malacca  and  Singapore  are  divided  in  to  nine  sectors  by 
Singapore’s  Vessel  Traffic  Information  System  (VTIS),  operated  by  Maritime  Port 
Authority  (MPA).  Figure  2  shows  the  VTIS  Operational  Area  Sectors  1  to  9,  starting 
from  Port  Klang  to  Singapore  Strait  leading  to  the  South  China  Sea.  Figure  3  shows  a 
close-up  of  Sectors  7  to  9  in  the  Singapore  Strait. 

With  its  close  proximity  and  accessibility  to  the  Strait  of  Singapore,  Jurong  Island 
is  a  high-value  target  to  terrorists,  because  a  successful  attack  on  Jurong  Island  would 
cause  major  economic  and  energy  disruption  to  the  world  economy. 


9 


Figure  2.  VTIS  Operational  Area  Sectors  1  to  9  (From  [13]) 
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Figure  3.  VTIS  Operational  Area  Sector  7  to  9  (From  [13]) 
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C.  MARITIME  THREAT  CHARACTERISTICS 

The  information  in  this  section  comes  from  CRS  Report  for  Congress  [2].  The 
threat  of  maritime  terrorism  is  significant  and  varies  with  attack  scenarios.  These 
scenarios  can  be  categorized  according  to  five  common  dimensions:  perpetrators, 
objectives,  locations,  targets,  and  tactics.  These  dimensions  are  useful  for  discussing 
historical  instances  of  maritime  terrorism  and  for  defining  potential  scenarios  for  future 
maritime  attacks.  Table  1  can  aid  in  generating  numerous  scenarios  based  on  different 
combinations  of  the  dimensions. 

1.  Perpetrators 

It  is  important  to  identify  potential  perpetrators  for  evaluating  maritime  attacks  as 
the  perpetrators’  backgrounds  and  capabilities  will  bear  on  the  types  of  attacks.  For 
example,  terrorists  who  are  trained  in  seamanship  may  use  this  skill  to  ram  a  SAW  to  its 
target.  In  the  SAW  case,  the  perpetrators  are  thus  the  terrorists. 

2.  Terrorists’  Objectives 

The  objectives  of  maritime  terrorism  can  vary — from  causing  human  casualties  to 
economic  losses  or  other  negative  impacts.  The  consequences  of  a  maritime  terrorism 
attack  can  be  minor  or  major.  If  human  casualties  are  the  primary  objective,  the  potential 
targets  are  likely  to  be  passenger/cruise  ships.  If  economic  loss  is  the  principal  objective, 
terrorists  will  likely  pick  containers  ships,  oil  tankers,  or  ship  channels  as  targets  so  as  to 
disrupt  shipping  or  trading.  The  attack  on  Limburg  was  of  such  an  objective;  it  caused  a 
reduction  in  Yemeni  shipping  volume  by  50%  and  increased  shipping  insurance  premium 
by  300%.  Also,  an  undesirable  impact  of  the  Limburg  attack  was  the  spilling  of  90,000 
barrels  of  oil,  causing  environmental  damage  [1], 

3.  Locations  of  Attacks 

Terrorist  attacks  on  vessels  are  likely  to  occur  in  areas  of  high  shipping  activities, 
especially  in  ports  or  at  shipping  chokepoints,  such  as  the  Straits  of  Malacca  and 
Singapore. 
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4. 


Terrorists’  Targets 


A  potential  target  for  a  SAW  can  be  an  oil  tanker,  a  port,  or  a  shore -based 
petrochemical  plant. 

5.  Terrorists’  Tactics 

In  a  maritime  terrorism  attack,  terrorists  can  employ  simple  and  yet  effective 
tactics  such  as  ramming  a  maritime  target  with  an  explosives-laden  dinghy.  As 
aforementioned,  the  USS  Cole  and  Limburg  attacks  demonstrated  the  effectiveness  of 
such  a  tactic.  A  more  complex  tactic  is  to  hijack  a  large  vessel,  in  particular  a  vessel  with 
petroleum/chemical  cargo,  and  ram  it  to  another  high-value  vessel  or  a  critical  shore- 
based  target;  the  petroleum/chemical  cargo  serves  as  the  explosive  component  of  the 
attack. 


Table  1.  Example  of  maritime  threat  characteristics  (From  [2]) 


Dimensions 

Example  Characteristics 

Perpetrators 

•  A1  Qaeda  and  affiliates 

•  Islamist  unaffiliated 

•  Foreign  nationalists 

•  Disgruntled  employees 

•  Others 

Objectives 

•  Mass  casualties 

•  Port  disruption 

•  Trade  disruption 

•  Environmental  damage 

Locations 

•  Shipping  chokepoints 

•  Ports 

•  Shore  based  petrochemical  processing  infrastructure 

Targets 

•  Military  vessels 

•  Cargo  vessels 

•  Fuel  tankers 

•  Ferries  /  cruise  ships 

•  Port  area  populations 
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Dimensions 

Example  Characteristics 

•  Ship  channels 

•  Port  industrial  plants 

•  Offshore  platforms 

Tactics 

•  Explosives  in  suicide  boats 

•  Explosives  in  light  aircraft 

•  Ramming  with  vessels 

•  Ship-launched  missiles 

•  Harbor  mines 

•  Underwater  swimmers 

•  Unmanned  submarine  bombs 

•  Exploding  fuel  tankers 

•  Explosives  in  cargo  ships 

•  WMDs  in  cargo  ships 

D.  SCENARIO  DEFINITION 

The  scenario  defined  in  this  research  is  now  described.  A  group  of  terrorists, 
whose  objective  is  to  disrupt  the  global  oil  trade  by  destroying  one  of  the  shore-based 
petrochemical  processing  plants  located  at  the  southernmost  end  of  Jurong  Island,  hijacks 
Ocean  Jewel  [15],  an  oil  tanker  on  its  westbound  shipping  course  in  Sector  7  of  the 
Singapore  Strait  (Figure  5).  As  the  hijacked  vessel  reaches  the  point  at  the  shortest 
distance  to  the  targeted  plant,  the  terrorists  maneuver  it  into  an  abrupt  right  turn  to  head 
towards  the  plant  at  full  speed  with  the  intent  to  ram  it  to  the  plant,  exploiting  its  cargo  of 
petroleum  product  as  the  explosive  component  in  order  to  magnify  the  damage  that  the 
attack  can  inflict.  The  oil  tanker  in  the  hands  of  the  terrorists  has  thus  become  a  SAW. 

Based  on  Table  2,  which  shows  both  the  SAW  (Ocean  Jewel)  details  and  the 
operating  conditions  for  the  scenario,  the  SAW  attack  begins  at  4  km  from  its  intended 
target,  and  the  estimated  time  of  SAW  collision  with  its  intended  target  is  5  minutes.  The 
ship  tenninology  used  in  Table  2  is  explained  in  Figure  4.  The  length  overall  (LOA) 
refers  to  the  distance  between  the  ship  extremities.  The  length  between  perpendiculars 
(LBP)  is  the  distance  between  the  after  and  forward  perpendiculars,  measured  parallel  to 

the  load  line  [14].  The  length  at  the  waterline  (LWL)  is  the  distance  from  the  forward 
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most  point  of  the  waterline  measured  in  profile  to  the  stem-most  point  of  the  waterline. 
The  beam  is  the  breadth  of  a  ship  at  the  widest  point.  The  depth  is  the  vertical  distance 
from  the  lowest  point  of  the  hull  to  the  deck  level.  The  draft  is  the  vertical  distance  from 
the  lowest  point  of  the  hull  to  the  water  level.  The  freeboard  is  the  difference  between 
the  depth  and  the  draft. 


Table  2.  Summary  of  SAW  details  (From  [15])  and  operating  conditions  (From  [7]) 


SAW  Characteristics: 

Vessel  Type 

Crude  oil  tanker 

Length  between  perpendiculars  - 
LBP  (m) 

265 

Beam  (m) 

43.2 

Depth  (m) 

23.8 

Draft  (m) 

17.38 

Full  Speed  (kts) 

25 

Full  Speed  (m/s) 

12.86 

Operating  Conditions: 

Sea  State 

<2 

(calm  rippled 
sea  with  a  wave 
height  no  more 
than  0.1  meters) 

Density  of  Seawater  (kg/m3) 

1025 

Gravitational  acceleration  (m/s2) 

9.81 
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Figure  5.  The  East-to-West  route  taken  by  the  terrorists  to  attack  Jurong  Island  (From 

[13]) 


In  regards  to  the  counter-SAW  operational  aspect  of  the  scenario,  the  following 
assumptions  are  made: 

•  A  WPB  is  applied  only  once. 

•  The  SAW  does  not  maneuver  to  counter  the  effects  of  the  WPB 
application. 

•  The  duration  of  an  effective  WPB  application  is  as  small  as  possible  so  as 
to  minimize  the  burden  on  WPB  generation. 

The  SAW’s  destructive  course  must  be  altered  as  early  as  possible  and  a  quickly 
as  possible  after  the  WPB  application. 
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E. 


MISSION  ANALYSIS  RESULTS 


There  are  various  ways  to  neutralize  the  SAW  in  the  above  SAW  scenario.  But, 
as  aforementioned,  an  operational  requirement  is  that  a  measure  to  counter  the  SAW  must 
not  sink  or  destroy  the  SAW.  A  measure  that  satisfies  this  requirement  is  to  divert  the 
SAW  from  its  destructive  course  toward  its  targeted  plant  on  Jurong  Island. 

In  Figure  6,  the  red  line  indicates  the  SAW  destructive  course  toward  its  targeted 
plant.  The  green  lines  indicate  the  off-collision  course  required  to  steer  clear  of  the  target 
to  either  port  (left)  or  starboard  (right)  side  of  the  SAW.  With  the  standoff  distance  of 
four  kilometers,  a  deflection  angle  of  30°  is  thus  required. 


Figure  6.  Illustration  of  the  SAW  scenario  analysis  (From  [13]) 
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III.  WATER  PLUME  BARRIER 


This  chapter  discusses  the  formation  of  a  water  plume  barrier  (WPB)  and  its 
potential  application  in  a  counter-SAW  mission. 

A.  WATER  PLUME  BARRIER  FORMATION 

A  WPB  is  the  free  water  surface  that  is  being  pushed  up  by  the  bubbles  created 
by,  for  example,  underwater  explosions.  A  set  of  bubbles  initiated  under  high  pressure 
will  oscillate  and  eventually  collapse  asymmetrically  with  the  water  jet  directed  away 
from  the  initial  quiescent  free  surface  [17].  As  a  result,  a  WPB  is  fonned. 

It  is  of  interest  to  create  in  a  short  time  a  WPB  shape  with  a  determined  series  of 
distributed  bubbles  created  by  underwater  explosions.  Certain  functionality  requires 
certain  shape  of  the  water  plume  barrier.  A  desired  shape  can  be  obtained  with  the  right 
values  of  the  bubble  parameters,  such  as  lateral  positions,  depths,  and  strengths.  To 
simplify  the  optimization  of  the  bubble  parameters,  the  Proper  Orthogonal 
Decomposition  (POD)  is  used  [17].  In  this  approach,  a  set  of  solution  snapshots  of  the 
bubble  and  free  surface  interaction  problem  are  used  to  fonn  a  set  of  linear  POD  basis 
functions,  which,  together  with  their  POD  coefficients,  are  then  used  to  optimize  the 
bubble  parameters  in  order  to  attain  the  desired  WPB  shape.  A  detailed  discussion  of  the 
use  of  the  POD  method  to  optimize  the  bubble  parameters  can  be  found  in  [17]. 

As  an  example,  a  desired  3-D  WPB  is  shown  in  Figure  7.  The  3-D  WPB 
pertaining  to  10  distributed  underwater  explosion  bubbles,  shown  in  Figure  8,  is  created 
using  the  POD  approach.  This  approach  is  efficient;  achieving  this  3-D  WPB  takes  50 
seconds  on  a  Intel  Xeon  2.8GHz  processor,  RAM  2.0Gb  [17]. 

The  next  section  discusses  the  use  of  a  WPB  in  a  counter-SAW  mission. 
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Figure  7.  Desired  WPB  (From  [17]) 


Figure  8.  Constructed  WPB  using  ten  bubbles  (From  [  1 7]) 

B.  WPB  USED  IN  COUNTER-SAW  MISSION 

A  set  of  underwater  bubbles  created  by  underwater  explosions  can  lead  to  the 
formation  of  a  water  plume  in  the  sea  surface.  These  bubbles  are  strategically  located  to 
create  a  particular  shape  WPB  to  alter  the  SAW  destructive  course  in  the  yaw  direction. 
Figure  9  illustrates  a  set  of  underwater  bubbles  which  forms  the  WPB  and  the  interaction 
between  the  WPB  and  the  SAW  hull. 
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The  WPB  induces  a  force  on  the  SAW  hull,  and  this  WPB-induced  force  needs  to 
be  of  a  certain  minimum  magnitude  and  the  duration  of  its  application  needs  to  be 
sufficient  long  in  order  to  steer  the  SAW  clear  off  its  destructive  course  (a  deflection 
angle  of  30°  as  stated  in  Section  E  of  Chapter  II).  Chapter  VI  elaborates  the  minimum 
magnitude  of  the  force  required  and  the  duration  and  location  of  application  of  the  force. 
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IV.  PARAMETRIC  STUDY 


As  discussed  in  Chapter  I,  three  parameters  to  be  detennined  in  countering  a 
SAW  are:  the  magnitude  of  the  WPB-induced  force,  the  location  of  its  application  along 
the  hull  of  the  SAW,  and  the  duration  of  its  application.  As  discussed  in  Chapter  II,  a 
successful  counter-SAW  mission  requires  that,  at  a  standoff  distance  of  four  kilometers 
from  its  targeted  plant  on  Jurong  Island,  the  SAW  be  steered  by  30°  to  either  its  port 
(left)  or  starboard  (right)  side  off  its  destructive  course  toward  its  target.  The  objective  of 
the  parametric  study  conducted  in  this  research  is  to  detennine  a  combination  or 
combinations  of  these  three  parameters  that  satisfy  the  counter-SAW  mission  success 
requirement. 

This  chapter  captures  the  parametric  study.  It  involves  generating  the 
mathematical  problem  of  detennining  the  response  (including  rotation  and  translation)  of 
a  ship  exerted  by  a  WPB-induced  force  and  solving  the  resulting  mathematical  problem, 
using  MATLAB  to  obtain  the  SAW  motion  and  the  angles  of  deflection  corresponding  to 
various  combinations  of  the  three  parameters.  The  fonnulation  of  the  mathematical 
problem  of  determining  the  response  involves  setting  up  the  equations  of  motion 
governing  the  SAW  dynamics,  defining  initial  conditions,  detennining  the  parameters  in 
the  equations  of  motion,  such  as  the  ship  characteristics,  the  magnitude  of  the  applied 
WPB-induced  force,  and  the  ship  hydrodynamics  derivatives.  A  computer  program  is 
used  to  generate  the  ship  hydrodynamics  derivatives,  given  the  SAW  characteristics  and 
its  operating  environment.  Various  magnitudes  of  the  applied  WPB-induced  forces  are 
derived  from  a  nominal  WPB-induced  force,  which  is  estimated  using  ship 
hydrodynamics. 

The  input  to  this  parametric  study  is  tied  to  the  parameters  in  the  equations  of 
motion.  These  parameters  are  divided  into  two  groups:  fixed  parameters  and  variable 
parameters.  The  fixed  parameters  are  the  ship  characteristics  and  the  ship  hydrodynamics 
derivatives.  The  variable  parameters  are  the  magnitudes  of  the  applied  WPB-induced 
forces,  the  durations  of  their  applications,  and  the  locations  of  their  applications  along  the 
hull  of  the  SAW. 
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The  output  of  this  parametric  study  consists  of  the  temporal  evolution  of  the  SAW 
translational  and  rotational  variables  and  the  deflection  angle  of  the  SAW  path  for 
various  combinations  of  the  magnitude  of  the  WPB-induced  force,  the  location  of  its 
application  along  the  hull  of  the  SAW,  and  the  duration  of  its  application  of  the  WPB- 
induced  force. 

A.  MATHEMATICAL  FORMULATION  OF  SAW  RESPONSE 
1.  Equations  of  Motion 

Two  coordinate  systems  are  used:  the  earth-fixed  coordinate  system  and  the  ship- 
fixed  coordinate  system  (Figure  10).  The  xo-axis,  yo-axis,  and  zo-axis  of  the  earth-fixed 
coordinate  system  form  a  right-handed  coordinate  system.  The  x-axis,  y-axis,  and  z-axis 
of  the  ship-fixed  coordinate  system,  whose  origin  is  at  the  amidships  (AMS),  form  a 
right-handed  coordinate  system. 


Figure  10.  Earth- fixed  coordinate  system  and  the  ship-fixed  coordinate  system. 


Next,  the  equations  of  motion  in  the  horizontal  plane  (i.e.,  the  xy-plane)  that 
govern  the  motion  of  a  SAW  (or  any  ship),  which  are  derived  in  [18],  are  captured  in  (1). 

(m  -  Y.  )v  -  (Y.  -  mxG)r  =  Yvv  +  (V  -  mU)r  +  FWPB  +  YsS 
(/.  -  N,  )r  -  (A,  -  mxG  )v  =  Nvv  +  (Nr  -  mxGU)r  +  NWPB  +  NsS 
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In(l), 


M :  the  vessel’s  mass 

v  :  the  velocity  in  the  x-direction  (known  as  the  sway  velocity) 

r  :  the  angular  velocity  in  the  clockwise  direction  (known  as  the  yaw  angular 
velocity) 

v :  the  acceleration  in  the  y-direction 
r  :  the  angular  acceleration  in  the  clockwise  direction 
7Z :  the  moment  of  inertia  about  the  z-axis 

xg  ■  the  distance  along  the  x-axis  from  the  amidships  (AMS)  to  the  ship’s  center 
of  gravity  (CG) 

Fwpb:  the  magnitude  of  the  WPB-induced  force 

Nwpb:  the  magnitude  of  the  torque  due  to  the  WPB-induced  force 

U :  the  nominal  forward  speed 

y :  the  change  in  force  in  the  sway  direction  per  a  unit  change  in  the  sway 
acceleration 

7  :  the  change  in  force  in  the  sway  direction  per  a  unit  change  in  the  sway 
velocity 

Y. :  the  change  in  force  in  the  sway  direction  per  a  unit  change  in  the  yaw  angular 
acceleration 

7  :  the  change  in  force  in  the  sway  direction  per  a  unit  change  in  the  yaw  angular 
velocity 

Ys :  the  change  in  force  in  the  sway  direction  per  a  unit  change  in  the  rudder 
deflection  angle 

Nv :  the  change  in  moment  in  the  yaw  rotation  per  a  unit  change  in  the  sway 
acceleration 
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Nv:  the  change  in  moment  in  the  yaw  rotation  per  a  unit  change  in  the  sway 
velocity 

N- :  the  change  in  moment  in  the  yaw  rotation  per  a  unit  change  in  the  yaw 
angular  acceleration 

Nr :  the  change  in  moment  in  the  yaw  rotation  per  a  unit  change  in  the  yaw 
angular  velocity 

Ns :  the  change  in  moment  in  the  yaw  rotation  per  a  unit  change  in  the  rudder 
deflection  angle 

The  quantities  with  the  subscripts  v,  r,  v;  r  and  5  are  the  hydrodynamics 
derivatives.  The  equations  of  motion  (1)  are  now  turned  into  a  system  of  linear 
differential  equations  for  the  velocity  v  and  the  angular  velocity  r;  that  is,  r  and  v  are 
linear  functions  of  r  and  v. 

v.  _  (4  ~  NrK  +  (K  ~ mxG)Nv  y  |  (4  ~ N,)(Yr  - m)  +  (Y.  - mxG)(Nr  - mxG)  ^ 

(4  - N.)(m  -  Y.)  -  {Y.  - mxG)(Nv  - mxG)  (/.  -  N,)(m  -Y.)- (Y.  - mxG)(N.  -  mxG) 

,  (I,  - Nt)Ys  +  (Yf  - mxG)Ns  §  |  (/,  - N,)FWPB  +  (Y.  - mxG)NWPB 

(4  -  N, )(m  -Y.)-  {Y.  -  mxG )(N.  -  mxG )  (/.  -  N, )(m  -Yv)-(Yr  -  mxG )(N.  -  mxG ) 

(2) 

■  =  (TV,  -  mxG)Yv  +  (m  -  Yi,)Nv  ^  |  (N,  -  mxG)(Yr  -  m)  +  (in  -  4)(7V,.  -  mxG)  ^ 

(I:  - N-)(m  -  Y.)  -  ( T :  - mxG)(Nv  - mxG)  (/.  - N-)(m  -Y.)-  ( T :  - mxa)(N ■  - mxG) 

_l _ (^Vy  ~  (m  ~  4)^Vg _ ^_| _ (2fy,  ~  )Fwpb  +  (m  —  YY)Nwpb _ 

(4  -  N,  )(m  -Y)-(Y  -  mxG  ){N.  -  mxG )  (Iz  -  N,  ){m  -Y.)-(Y.  -  mxG  )(N.  -  mxG ) 

Note  that,  in  this  work,  no  environmental  disturbances,  such  as  wind,  currents,  or 
waves,  are  considered,  and  the  rudder  deflection  angle  is  zero  degree.  Using  the  water 
density  p,  the  ship  length  L,  and  the  nominal  forward  speed  U,  the  variables  can  be 
nondimensionalized  as  shown  in  (3).  The  dimensionless  fonn  of  the  equations  of  motion 
in  the  primed  variables  turns  out  to  be  identical  to  that  of  (1)  and  (2)  with  IT= 1.  Unless 
stated  otherwise,  this  work  assumes  the  dimensionless  form  of  the  equations  of  motion,  , 
and,  for  convenience,  the  primes  are  dropped. 
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The  evaluation  of  the  hydrodynamics  derivatives,  which  are  dimensionless,  is 
discussed  in  Section  A  of  Chapter  VI. 


2.  Parameters 

The  equations  of  motion  in  (1)  or  the  equations  in  (2)  contain  two  types  of 
parameters:  fixed  and  variable  parameters.  The  variable  parameters  are  the  magnitude  of 
FWpb  and  the  magnitude  of  the  torque  Nwpb,  whose  value  depends  on  the  location  of 
application  of  FWPB  and  the  magnitude  of  FWPB  (the  angle  of  application  being  90°). 
Different  values  of  the  magnitude  of  Fwpb  are  based  on  a  nominal  value,  whose 
estimation  using  ship  hydrostatics  is  elaborated  in  Section  B  of  Chapter  IV.  While  not  a 
parameter  inherent  in  the  equations  of  motion,  the  duration  of  application  of  FWPB  is  also 
considered  in  the  study  of  the  response  of  the  SAW  to  the  application  of  FWPB  over  a 
small,  finite  duration  of  time.  The  parametric  nature  of  the  study  thus  reflects  the 
variation  of  the  magnitude  of  FWpb,  the  location  of  application  of  Fwpb,  and  the  duration 
of  application  of  Fwpb- 

The  fixed  parameters  are  the  hydrodynamics  derivatives,  the  SAW  mass,  m,  the 
moment  of  inertia  about  the  z-axis,  7Zj  and  the  distance  along  the  x-axis  from  the 
amidships  (AMS)  to  the  ship’s  center  of  gravity  (CG),  xg-  The  latter  three  parameters  are 
related  to  the  SAW  characteristics.  The  hydrodynamics  derivatives  are  computed  using 
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the  maneuvering  prediction  program  (MPP1.3)  [19],  which  uses  as  input  the  vessel  and 
operating  environment  characteristics  described  in  Chapter  II.  Appendix  A  contains  the 
details  of  MPP1.3. 

3.  Ship  Kinematics 

Given  an  initial  condition  of  r  and  v-  integration  of  the  system  in  (2)  yields 
r,v,r  and  v  as  functions  of  time  t.  The  yaw  angle,  y/ ,  is  measured  counterclockwise 
with  respect  to  the  xo-axis.  It  can  be  obtained  by  integrating  (//  ,  the  yaw  velocity  (in  the 
clockwise  direction),  with  respect  to  time  t.  In  particular,  over  a  sufficiently  small  time 
interval  [t,t  +  At],  y/  is  given  by 

y/{t  +  At)-y/(t)  +  y/{t)At  (4) 

or,  since  y/  is  equal  to  r, 

y/(t  +  At)-y/(t)  +  r(t)At  (5) 

Finally,  the  kinematical  variables  in  the  earth- fixed  frame,  xo,yo,x0  and  y0,  are 
obtained  according  to 

xQ  ( t )  =  u{t)  cos  y/(t)  -  v{t)  sin  y/{t) 
ya  (0  =  «(0  sin  y/{t)  +  v(t)  cos  y/(t) 

(6) 

x0(t  +  At)  =  x0(t)  +  x0(t)At 
yo  (t  +  At)  =  yo  (t)  +  yo  (t  )At 

where  u  is  the  velocity  in  the  y-direction  (known  as  the  surge  velocity). 

B.  NOMINAL  WPB-INDUCED  FORCE  CALCULATION 

A  simplified  representation  of  a  stationary  ship  in  water  is  a  rectangular  block 
immersed  in  undisturbed  water,  assumed  to  be  incompressible,  non-viscous,  and  with  no 
surface  tension.  The  ship  is  in  an  equilibrium  state  when  the  drafts  at  forward,  aft,  port 
and  starboard  are  equal. 

The  free-body  diagram  in  Figure  1 1  shows  the  hydrostatic  forces  exerted  on  the 
ship  in  its  equilibrium  state.  The  weight  of  the  ship  is  denoted  by  W.  The  buoyancy  force 
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(F Buoyancy)  's  the  upward  vertical  force  exerted  by  the  water;  its  magnitude  equals  to  the 

weight  of  the  fluid  displaced.  The  horizontal  forces  ( FS3  and  FS4 )  exerted  by  the  water 

are  of  equal  magnitudes  but  point  in  the  opposite  direction.  The  net  force  in  the 
horizontal  direction  is  thus  zero.  The  magnitudes  of  these  forces  are  calculated  according 
to  (7). 

FS3  =  FS4  =  L  Jo  pgzdz  =  ^ LpgT 2 
^Buoyancy  =  pgLBT 

where  L  is  the  length  of  the  ship,  B  the  breath  of  the  ship,  T  the  draft  of  the  ship,  z 
distance  from  the  water  surface  to  the  center  of  pressure  (or  the  depth  of  pressure),  p 
density  of  the  seawater,  and  g  the  gravitational  acceleration. 


(V) 

the 

the 


Figure  1 1 .  Hydrostatic  force  on  a  ship  body  at  equilibrium 


Figure  12  depicts  the  changes  in  the  hydrostatic  forces  when  a  WPB  is  applied  to 
the  ship.  As  the  WPB  is  applied,  the  forward  draft  (I)  at  the  starboard  would  be 

increased  by  aT,  for  0  <  a  <  1,  thereby  creating  an  additional  force,  F  ,  in  the  negative  y- 
direction.  The  resultant  horizontal  force,  P'S3 ,  which  is  the  sum  of  PS3  and  P  ,  is  now 
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greater  than  FS4 ,  causing  the  ship  to  roll  by  an  angle  O  and  resulting  in  FlB  ,  which 
deviates  from  the  z-axis  by  the  same  angle,  O. 


Fy 


Figure  12.  Changes  in  the  hydrostatic  forces 


With  the  application  of  the  WPB, 
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The  WPB  is  assumed  to  be  applied  at  some  location  in  the  forward  area  of  the 
ship.  The  footprint  of  the  WPB  on  the  SAW  (i.e.,  the  length  of  the  WPB  that  actually 
interacts  with  the  SAW)  covers  a  fraction  of  the  hull  of  the  SAW,  denoted  by  /?,  which 
takes  values  in  (0,  1).  The  magnitude  of  the  resultant  hydrostatic  force,  or  the  WPB- 

induced  force,  VWPK ,  exerted  on  the  hull  is  then  given  by 


F wpb  =  Fy/?  =  F\uoyancy,  sin  0/7  =  i  pgLB[T  +  T(  1  +  a)]  sin  O/i  (9) 
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This  VWPH  is  the  nominal  WPB-induced  force  to  be  used  in  the  generation  of 
variants  of  the  WPB-induced  forces  considered  in  the  parametric  study. 

The  location  of  application  of  the  WPB-induced  force,  assumed  to  be  in  the  fore 

section  of  the  hull,  is  at  a  distance  of  yL  from  the  bow,  where  0  <  y  <  — ,  and  a  distance  / 

(the  moment  arm  of  FWfVi )  from  the  amidships  (AMS).  As  the  geometry  in  Figure  13 


shows,  / 


I  > 

—  Y 

v2 


L  .  The  magnitude  Nwpb  of  the  moment  resulting  from  the  application 


of  the  WPB-induced  force,  FWPB ,  perpendicular  to  the  hull  is  then  given  by 


N  =IF 

1  v  WPB  11  WPB 


—  Y 

2  j 


LR 


WPB 


(10) 


Figure  13.  Plan  view  of  the  SAW  and  location  of  the  applied  FWPB 
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c. 


SOLUTION  IMPLEMENTATION 


1.  Calculation  of  Parameters 

Tables  3  and  4  summarize  the  input  parameters.  Two  counter-SAW  operational 
assumptions  are  made  in  the  determination  of  the  magnitude  of  the  nominal  WPB- 
induced  force. 

1)  The  factor  ot  (introduced  in  Chapter  IV)  by  which  the  SAW  draft  is  increased 
by  the  application  of  a  WPB  is  set  to  be  0.3. 

2)  The  footprint  p  is  set  to  be  0.3. 

With  the  roll  angle  <D>  obtained  from  (7)  and  the  values  of  p,  g,  L,  B,  and  T  given 
in  Table  2  in  Chapter  II,  the  calculation  of  the  magnitude  of  the  nominal  WPB-induced 
force,  Fwpb  ,  using  (8),  results  in  FWpb  =  82.7  x  106  N. 

The  location  of  VWPB  application  varies  between  0.1L  and  0.25L  from  the  bow  tip. 
The  duration  of  VWPB  application  varies  between  1  to  2.5  seconds. 

The  values  of  the  hydrodynamics  derivatives  in  Table  4  are  generated  by  the 
maneuvering  prediction  program  (MPP1.3)  mentioned  in  Section  A  of  Chapter  IV.  The 
values  of  the  derivatives  are  unitless. 


Table  3.  Values  for  Fwpb,  location  and  duration  of  Fwpb  application 


S/No. 

INPUT  PARAMETERS 

NOTATION 

VALUE 

1. 

Magnitude  of  WPB-induced  force  (N) 

Fwpb 

82.7  x  106 

to  206.8  x  106 

2. 

Location  of  Fwpb  application  (m) 

/ 

O.lLto  0.25L 

3. 

Duration  of  Fwpb  application  (s) 

1  to  2.5 

4. 

Dimensionless  mass 

m 

0.017620 

5. 

Dimensionless  moment  of  inertia 

Iz 

0.001101 
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Table  4.  Values  of  hydrodynamic  derivatives 


S/No. 

INPUT  PARAMETERS 

NOTATION 

VALUE 

1. 

The  change  in  force  in  the  sway  direction 
per  a  unit  change  in  the  sway  velocity 

V 

-0.025566 

2. 

The  change  in  force  in  the  sway  direction 
per  a  unit  change  in  the  sway  acceleration 

n 

-0.016110 

3. 

The  change  in  force  in  the  sway  direction 
per  a  unit  change  in  the  yaw  angular 
velocity 

K 

-0.008402 

4. 

The  change  in  force  in  the  sway  direction 
per  a  unit  change  in  the  yaw  angular 
acceleration 

Yr 

-0.001046 

5. 

The  change  in  moment  in  the  yaw  rotation 
per  a  unit  change  in  the  sway  velocity 

K 

0.005079 

6. 

The  change  in  moment  in  the  yaw  rotation 
per  a  unit  change  in  the  sway  acceleration 

N, 

-0.001200 

7. 

The  change  in  moment  in  the  yaw  rotation 
per  a  unit  change  in  the  yaw  angular 
velocity 

K 

-0.003691 

8. 

The  change  in  moment  in  the  yaw  rotation 
per  a  unit  change  in  the  yaw  angular 
acceleration 

N, 

-0.000870 

9. 

The  change  in  force  in  the  sway  direction 
per  a  unit  change  in  the  rudder  deflection 
angle 

Ys 

0.003275 

10. 

The  change  in  moment  in  the  yaw  rotation 
per  a  unit  change  in  the  rudder  deflection 
angle 

Ns 

-0.001605 

2.  MATLAB  Implementation 

The  equations  of  motion  in  (2),  (5)  and  (6)  are  coded  using  MATAB,  and  the 
resulting  code  is  run  with  the  values  of  the  parameters  in  Tables  3  and  4  to  generate  the 
results  for  the  parametric  analysis.  Appendix  B  contains  the  MATLAB  code. 
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a. 


MATLAB  Inputs 


(1)  Code  Execution  Data.  The  code  is  executed  for  five  seconds  with  a 
time  step  At  of  0.001  seconds. 

(2)  Hydrodynamic  Derivatives.  Table  4  contains  the  values  of  the 
hydrodynamics  derivatives  obtained  with  MPP1.3,  as  discussed  in  Section  A  of  Chapter 
IV.  The  vessel  and  operating  environment  characteristics  described  in  Chapter  II  are 
inputs  to  MPP1.3. 

(3)  Nominal  FWpb-  As  shown  in  Table  3,  the  nominal  magnitude  of  the 
WPB-induced  force  FWPB  is  82.7  x  106N. 

(4)  Initial  Conditions.  The  initial  conditions  refer  to  the  kinematical 
variables  at  the  time  immediately  prior  to  the  application  of  Fwpb- 

(5)  Parametric  Study  Data.  Four  different  magnitudes  of  the  WPB- 
induced  force,  Fwpb,  are  considered:  82.7  x  106  N,  corresponding  to  the  nominal  value, 
82.7  x  106  N  (Force  #1),  124.1  x  106  N  (Force  #2),  165.4  x  106  N  (Force  #3),  and  206.8  x 
106  N  (Force  #4).  Each  of  these  forces  can  be  applied  to  four  different  locations  along 
the  hull;  they  correspond  to  four  different  values  of  y:  0.1,  0.15,  0.2,  and  0.25.  For  each 
of  these  forces  and  each  location  of  its  application,  four  different  durations  of  the  WPB- 
induced  force  application  are  considered:  1,  1.5,  2,  and  2.5  seconds.  The  application  of 
these  forces  always  starts  at  1.0  second. 

b.  MATLAB  Outputs 

For  each  combination  of  the  three  parameters  —  the  magnitude  of  the 
WPB-induced  force,  the  location  of  its  application  along  the  hull  of  the  SAW,  and  the 
duration  of  its  application  of  the  WPB-induced  force,  the  following  outputs  are  obtained 
and  displayed  graphically. 

(1)  Angular  Deviation  of  SAW  Path.  The  angle  9  between  the  SAW’s 
intended  destructive  course  and  its  deflected  path  resulting  from  the  application  of  WPB 
is  computed  as  a  function  of  time.  The  deflected  path  is  along  the  bearing  direction  when 
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the  application  of  FWPB  stops.  This  result  is  used  in  the  assessment  of  the  counter-SAW 
mission  success. 

(2)  SAW  Motion  in  Ship-fixed  Coordinate  System.  The  SAW 

acceleration  v  in  the  sway  direction,  the  angular  acceleration  r  in  the  clockwise 
direction,  the  velocity  v  in  the  y-direction  are  obtained  as  functions  of  time.  These 
outputs  are  to  illuminate  the  SAW  motion  in  both  coordinate  systems. 
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V.  PARAMETRIC  STUDY  RESULTS  AND  DISCUSSION 


As  previously  discussed,  the  three  parameters  that  are  varied  in  this  parametric 
study  are:  the  magnitude  of  the  WPB-induced  force,  the  location  of  its  application  along 
the  hull  of  the  SAW,  and  the  duration  of  the  application  of  the  WPB-induced  force.  The 
objective  of  the  parametric  study  is  to  determine  a  combination  or  combinations  of  these 
three  parameters  that  satisfy  the  counter-SAW  mission  success  requirement. 

As  discussed  in  Chapter  IV,  r,r(t),v,v(t )  and  9  stand  for  the  angular  acceleration, 
the  angular  velocity,  the  acceleration  in  the  sway  direction,  the  velocity  in  the  y-direction, 
and  the  deflection  angle,  respectively.  The  parametric  study  results  of  { r,r(t),v,v(t )  and 
0}  are  obtained  for  the  following  three  combinations  of  the  parameters. 

1)  Variable  FWpb,  a  fixed  location  of  FWPB  application,  and  a  fixed  duration  of  FWPB 

application. 

2)  A  fixed  Fwpb,  variable  locations  of  ¥WPB  application,  and  a  fixed  duration  of 

Fwpb  application. 

3)  A  fixed  Fwpb1  a  fixed  location  of  FWPB  application,  and  variable  durations  of 

VWPB  application. 

A.  COMBINATION  1— VARIABLE  MAGNITUDE  OF  %PB 

Figures  14  and  15  show  the  temporal  evolution  of  v  and  r ,  respectively,  for  the 
four  indicated  forces,  the  location  and  duration  of  ¥WPB  application  remaining  fixed  at  y  = 
0.15  and  one  second,  respectively.  The  sway  acceleration,  v,  and  the  angular 
acceleration,  r ,  reach  their  peak  values  when  each  of  these  forces  is  applied  and,  because 
the  resistance  of  the  hydrodynamic  force  overcomes  FWPB,  they  rapidly  decrease  until  the 

F WPB  application  ends,  at  which  time  they  peak  again  but  in  the  opposite  direction.  They 
then  decrease  rapidly  to  zero  due  to  the  damping  effects.  The  damping  effects  refer  to  the 
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effects  of  the  resistance  of  the  hydrodynamic  force.  Furthennore,  the  largest  sway 
acceleration  and  the  largest  angular  acceleration  correspond  to  the  largest  magnitude  of 

VWPB  (i.e.,  Force  #4). 


Time  vs  vdot 


Figure  14.  Time  vs  v 


Time  vs  rdot 


Figure  15.  Time  vs  f 


36 


Figures  16  and  17  show  the  temporal  evolution  of  v  and  r,  respectively,  for  the 
four  indicated  forces,  the  location  and  duration  of  VWPB  application  remaining  fixed  at  y  = 

0.15  and  one  second,  respectively.  During  the  YWPB  application,  the  velocity,  v,  and  the 

angular  velocity,  r,  rapidly  reach  their  peak  values  and,  when  the  VWPB  application  ends, 
they  decrease  rapidly  to  zero  because  of  the  damping  effects.  Furthermore,  the  largest 
velocity  and  angular  velocity  correspond  to  the  largest  magnitude  of  VWPB  (i.e.,  Force  #4). 


Time  vs  v 


Figure  16.  Time  vs  v 
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Time  vs  r 


Figure  17.  Time  vs  r 

Figure  18  shows  the  temporal  evolution  of  the  deflection  angle  9  for  the  four 
indicated  forces,  the  location  and  duration  of  Fwpb  application  remaining  fixed  at  y  =  0. 15 
and  one  second,  respectively.  It  takes  2.3  seconds  to  achieve  the  required  deflection 
angle  of  30°  with  the  largest  magnitude  of  the  WPB-induced  force  and  about  4.5  seconds 
with  its  smallest  magnitude. 
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Angle  of  deflection  vs  time 


Figure  18.  Angle  of  deflection  vs  time 

This  set  of  results  thus  indicates  that  the  high  magnitude  of  the  WPB-induced 
force  is  more  effective  in  diverting  the  destructive  course  than  does  to  the  low  magnitude. 
The  required  deflection  angle  of  30°  is  achieved  more  quickly  with  the  high  magnitude 
than  with  the  low  magnitude  of  the  WPB-induced  force. 

B.  COMBINATION  2— VARIABLE  LOCATION  OF  FWPB  APPLICATION 

Figures  19  shows  the  temporal  evolution  of  the  sway  acceleration,  v,  for  the  four 
indicated  locations  of  FWPB  application,  the  magnitude  of  the  applied  Fwpg  and  the 
duration  of  FWPB  application  remaining  fixed  at  82.7  x  106  N  and  one  second, 
respectively.  The  temporal  behavior  of  vis  identical  to  that  of  v  in  Figure  14, 
irrespective  of  the  location  of  FWPB  application,  in  the  sense  that  v  reaches  its  peak  value 

when  Fwpb  is  applied  and,  because  the  resistance  of  the  hydrodynamic  force  overcomes 

the  Fwpb,  it  rapidly  decreases  until  the  FWPB  application  ends,  at  which  time  it  peaks  again 
but  in  the  opposite  direction.  It  then  decreases  rapidly  to  zero  due  to  the  damping  effects. 
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As  the  location  of  YWPB  application  recedes  from  the  bow  (i.e.,  y  increases),  the  applied 

moment  NWpb  decreases.  The  largest  moment  results  in  the  smallest  sway  acceleration  as 
the  translation  effect  is  at  the  minimum.  However,  as  shown  in  Figure  20,  the  largest 
angular  acceleration,  r  ?  corresponds  to  the  largest  moment  as  it  produces  the  largest 
rotation  effect. 


Time  vs  vdot 


Figure  19. 


Time  vs  v 
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3 


Time  vs  rdot 


Figure  20.  Time  vs  r 


Figures  21  and  22  show  the  temporal  evolution  of  v  and  r,  respectively,  for  the 
four  indicated  locations  of  VWPB  application,  the  magnitude  of  the  applied  VWPB  and  the 

duration  of  fWPB  application  remaining  fixed  at  82.7  x  106  N  and  one  second, 
respectively.  During  the  F^  application,  the  velocity,  v,  and  the  angular  velocity,  r, 

rapidly  reach  their  peak  values  and,  when  the  VWPB  application  ends  (i.e.,  the  application 
of  the  moment  ends),  they  decrease  rapidly  to  some  values  because  of  the  damping 
effects.  As  the  location  of  VWPB  application  recedes  from  the  bow  (i.e.,  y  increases),  the 

applied  moment  NWpb  decreases.  The  largest  moment  results  in  the  smallest  velocity  as 
the  translation  effect  is  at  the  minimum.  However,  as  shown  in  Figure  22,  the  largest  the 
angular  velocity,  r,  corresponds  to  the  largest  moment  as  it  produces  the  largest  rotational 
effect.  The  smaller  the  moment  is,  the  greater  the  SAW  restoring  moment  is;  the 
resulting  effect  is  an  increase  of  the  angular  velocity  in  the  positive  region. 
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Time  vs  v 


Figure  2 1 .  Time  vs  v 


Time  vs  r 


Figure  22.  Time  vs  r 
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Figure  23  shows  the  temporal  evolution  of  the  deflection  angle  9  for  the  four 
indicated  locations  of  VWPB  application,  the  magnitude  of  the  applied  fWPB  and  the  duration 

of  ^wpb  application  remaining  fixed  at  82.7  x  106  N  and  one  second,  respectively.  The 
location  ofFW7Vi  application  closest  to  the  bow  (y  =  0.1),  hence  the  largest  moment,  allows 
the  required  deflection  angle  of  30°  to  be  achieved  in  3.6  seconds,  the  smallest  amount  of 
time  as  compared  to  those  corresponding  to  the  other  locations  of  VWPB  application. 


Angle  of  deflection  vs  time 


Figure  23.  Angle  of  deflection  vs  time 

This  set  of  results  thus  indicates  that  the  greater  moments  are  more  effective  in 
diverting  the  destructive  course  than  do  the  greater  moments.  The  required  deflection 
angle  of  30°  is  achieved  more  quickly  with  the  moment  arms  corresponding  to  y  =  0.1 
and  0. 15  than  with  the  other  moment  arms. 

C.  COMBINATION  3— VARIABLE  DURATION  OF  VWPB  APPLICATION 

Figures  24  and  25  show  the  temporal  evolution  of  v  and  r ,  respectively,  for  the 
four  indicated  durations  of  VWPB  application,  the  magnitude  of  Fj(7V;  and  the  location  of 
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¥wpb  application  remaining  fixed  at  82.7  x  106  N  and  y  =  0.15,  respectively.  The  sway 

acceleration,  v ,  and  the  angular  acceleration,  r ,  reach  their  peak  values  when  each  of 
these  forces  is  applied  and,  because  the  resistance  of  the  hydrodynamic  force  overcomes 
PWPB  ,  they  rapidly  decrease  until  the  FWPB  application  ends,  at  which  time  they  peak  again 

but  in  the  opposite  direction.  They  then  decrease  rapidly  to  zero  due  to  the  damping 
effects.  Furthermore,  the  times  at  which  they  peak  in  the  opposite  direction  naturally 

reflect  the  durations  of  FWPB  application.  The  temporal  evolution  profiles  of  v  and  r  are 
identical,  irrespective  of  the  durations  of  FWPB  application,  because  the  magnitude  of  FWPB 
and  the  location  of  FWPB  application  are  fixed. 


Time  vs  vdot 


Figure  24. 


Time  vs  v 
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2.5 


Time  vs  rdot 


Figure  25.  Time  vs  r 


Figures  26  and  27  show  the  temporal  evolution  of  v  and  r,  respectively,  for  the 
four  indicated  durations  of  VWPB  application,  the  magnitude  of  VWPB  and  the  location  of 

VWPR  application  remaining  fixed  at  82.7  x  106  N  and  y  =  0.15,  respectively.  During  the 

Fwpb  application,  the  velocity,  v>  and  the  angular  velocity,  r,  rapidly  reach  their  peak 
values  and,  when  the  FWpb  application  ends,  they  decrease  rapidly  to  zero  because  of  the 

damping  effects.  Furthermore,  the  larger  the  duration  of  VWPB  application,  the  longer  the 
duration  of  the  rapid  increase  of  v  and  r  during  the  acceleration  phase,  and  the  smaller  the 
time  it  takes  v  and  r  to  rapidly  decrease  during  the  deceleration  phase. 
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Time  vs  v 


Figure  26.  Time  vs  v 


Time  vs  r 


Figure  27.  Time  vs  r 
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Figure  28  shows  the  time  that  is  required  to  achieve  the  deflection  angle  of  30°, 
for  the  four  indicated  durations  of  VWPB  application,  the  magnitude  of  VWPB  and  the 

location  of  FWPB  application  remaining  fixed  at  82.7  x  106  N  and  y  =  0.15,  respectively. 

For  the  2.5-second  duration  of  FWPB  application,  it  takes  about  3.2  seconds  to  reach  the 
required  deflection  angle.  An  amount  of  4.5  seconds  is  needed  to  reach  required 
deflection  angle  for  the  1.0-second  duration  of  VWPB  application. 


Angle  of  deflection  vs  time 


Figure  28.  Angle  of  deflection  vs  time 


This  set  of  results  thus  indicates  that  the  long  durations  of  VWPB  application  are 

more  effective  in  diverting  the  SAW  destructive  course  than  do  the  short  durations.  The 
longer  durations  result  in  reaching  the  required  deflection  angle  faster. 

The  obtained  results  presented  in  this  chapter  are  used  in  the  preliminary 
assessment  of  the  feasibility  of  using  a  WPS  to  counter  a  SAW,  which  is  discussed  in  the 
next  chapter. 
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VI.  PRELIMINARY  ASSESSMENT 


An  assessment  of  the  feasibility  of  using  a  WPB  to  counter  a  SAW  amounts  to 
answering  to  the  questions  posed  in  the  problem  stated  in  Chapter  I:  Is  it  feasible  to 
apply  WPB-induced  forces  to  the  hull  of  the  SAW  to  alter  its  destructive  course?  If 
feasible,  what  are  the  magnitudes  of  the  WPB-induced  forces  and  for  how  long  do  they 
need  to  be  maintained?  Where  on  the  hull  of  the  SAW  should  they  be  applied?  Put 
differently,  which  combinations  of  the  three  parameters  —  the  magnitude  of  the  WPB- 
induced  force,  the  location  of  its  application  along  the  hull  of  the  SAW,  and  the  duration 
of  its  application  of  the  WPB-induced  force  —  satisfy  the  counter-SAW  mission 
requirement  of  steering  the  SAW  30°  off  its  destructive  course  established  in  Chapter  II. 
This  chapter  provides  this  assessment. 
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A. 


SUMMARY  OF  PARAMETRIC  STUDY  RESULTS 


Table  5.  Summary  of  parametric  study  results 


WPB-Induced  Force, 

Fwpb  (N) 

Location  of 

Fwpb 

Application,  y 

Duration  of 

Fwpb 

Application  (s) 

Time  to  Deflection 

Angle  of  30°  (s) 

Variable 

Fixed 

Fixed 

82.7  x  106 

0.15 

1 

4.5 

124.1  x  106 

3 

165.4  x  106 

2.5 

206.8  x  106 

2.3 

Fixed 

Variable 

Fixed 

82.7  x  106 

0.1 

1 

3.5 

0.15 

4.5 

0.2 

03 

0.25 

00 

Fixed 

Fixed 

Variable 

82.7  x  106 

0.15 

1 

4.5 

1.5 

3.4 

2 

3.2 

2.5 

3.2 

The  parametric  study  results  captured  in  Table  5  indicate  that  large  magnitudes  of 
F WPB  result  in  satisfying  the  counter-SAW  mission  success  requirement  of  30°  angle  of 

deflection  in  short  times.  A  VWPB  of  82.7  x  106  N  magnitude  applied  for  one  second  at  a 
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location  between  0.1  and  0.5  the  length  (L)  of  the  SAW  from  its  bow  takes  4.5  seconds  to 
achieve  the  deflection  angle  of  30°.  Since  the  SAW  moving  at  12.86  m/s  would  reach  its 
intended  target  on  Jurong  Island  in  at  least  five  minutes  from  a  standoff  distance  of  four 
kilometers  from  the  target,  this  combination  of  the  magnitude  of  the  WPB-induced  force, 
the  location  of  its  application  along  the  hull  of  the  SAW,  and  the  duration  of  its 
application  of  the  WPB-induced  force  should  sufficiently  satisfy  the  counter-SAW 
mission  success  requirement. 

This  combination  serves  as  the  preliminary  requirements  that  a  WPB  impart  a 
force  of  82.7  x  106  N  magnitude  that  is  exerted  for  one  second  at  a  location  along  the 
SAW  hull  somewhere  between  0.1L  and  0.1 5L  from  the  bow. 
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VII.  CONCLUSION  AND  RECOMMENDATIONS 


This  chapter  recapitulates  the  problem,  provides  a  summary  of  the  research 
findings,  and  recommendations  for  future  work. 

A.  CONCLUSION 

The  specific  SAW  problem  addressed  in  this  thesis  involves  an  oil  tanker,  Ocean 
Jewel,  commandeered  by  terrorists  to  ram  a  petrochemical  processing  plant  located  at  the 
southernmost  end  of  Jurong  Island  in  the  Singapore  Strait.  If  sunk  or  destroyed  as  a 
result  of  a  counter  measure,  the  SAW’s  collateral  damage  would  severely  disrupt  the 
traffic  flow  in  the  shipping  lane.  To  prevent  such  a  disruptive  catastrophe,  non¬ 
destructive  measures  must  be  implemented  to  cause  the  oil  tanker  to  deviate  from  its 
destructive  path  toward  the  target.  One  such  a  measure  involves  a  strategic  application  of 
WPB-induced  forces  to  the  SAW.  A  WPB  is  the  free  water  surface  that  is  being  pushed 
up  by  the  bubbles  created  by,  for  example,  underwater  explosions.  This  thesis  examines 
the  feasibility  of  using  a  strategically  created  WPB  to  alter  the  destructive  course  of  the 
SAW. 

The  feasibility  is  assessed  from  two  aspects:  (1)  The  effectiveness  of  a 
combination  of  the  magnitude  of  the  WPB-induced  force,  the  location  of  its  application 
along  the  hull  of  the  SAW,  and  the  duration  of  its  application  in  diverting  the  SAW  from 
its  path  toward  its  target  and  (2)  the  feasibility  of  generating  a  WPB  that  can  be  used  as 
specified  by  the  combination. 

The  research  captured  in  this  thesis  deals  mainly  with  the  first  feasibility  aspect. 
The  findings  in  this  research  indicate  that  it  is  feasible  to  use  a  WPB  to  counter  the  SAW. 
Specifically,  at  a  standoff  distance  of  four  kilometers  from  its  targeted  plant  on  Jurong 
Island,  the  SAW  can  be  diverted  by  30°  from  its  path  toward  the  target  by  a  WPB- 
induced  force  of  82.7  x  106  N  magnitude  exerted  for  one  second  at  a  location  along  the 
SAW  hull  somewhere  between  0.1L  and  0.15L  from  the  bow.  The  WPB  length  covers  a 
fraction  of  0.3  of  the  hull  of  the  SAW  and  the  force  increases  the  saw  draft  by  a  factor  of 
0.3. 
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The  feasibility  of  generating  a  WPB  that  can  produce  a  force  of  the  required 
magnitude  that  can  be  sustained  for  the  required  duration  is  currently  investigated  by 
National  University  of  Singapore. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

In  Section  C  of  Chapter  IV,  the  factor,  a,  by  which  the  SAW  draft  is  increased 
and  the  footprint  of  the  WPB,  P,  are  both  assumed  to  be  0.3.  A  future  effort  is 
recommended  to  determine  other  values  for  these  quantities  for  which  the  use  of  a  WPB 
to  counter  the  SAW  is  feasible.  Furthennore,  it  is  recommended  that  the  on-going  NUS 
research  detennine  the  feasibility  of  generating  water  plume  barriers  with  shapes  and 
induced  forces  that  satisfy  those  values  of  a  and  p. 

The  work  in  this  thesis  assumes  zero  rudder  deflection  angle,  no  environmental 
disturbances  such  as  wind,  current,  and  no  countermeasures  from  the  SAW  like  steering 
the  SAW  back  to  its  destructive  course.  A  modeling  and  simulation  effort  is 
recommended  to  take  into  account  non-zero  rudder  deflection  angle,  environmental 
disturbances,  and  countermeasures  from  the  SAW  in  the  assessment  of  the  counter-SAW 
mission  effectiveness.  The  results  of  this  effort  could  influence  system  concepts  and 
counter-SAW  operations. 
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APPENDIX  A.  DETERMINATION  OF  HYDRODYNAMIC 

DERIVATIVES 


A  maneuvering  prediction  program  (MPP1.3)  is  used  to  detennine  the 
hydrodynamic  derivatives.  This  Window-based  computer  program  is  developed  by  M.G. 
Parsons,  Department  of  Naval  Architecture  and  Marine  Engineering,  University  of 
Michigan  [19].  It  offers  two  main  options,  namely  the  linear  evaluation  for  the 
assessment  of  ship  course  stability,  turning  ability  and  controllability  and  the  turning 
prediction  for  the  estimation  of  turning  circle  characteristics.  In  this  thesis,  only  the  linear 
evaluation  is  used. 


1.  Input  to  Maneuvering  Prediction  Program  (MPP1.3) 

The  input  to  MPP  is  through  a  series  of  five  windows  within  the  menu 
inputs  as  follows: 

•  Project  Name 

•  Vessel  Characteristics 

•  Steering  Characteristics 

•  Operating  Conditions 

•  Water  Properties 

Project  Name  provides  a  location  for  the  general  identification  of  the 
project  being  analyzed. 


Figure  29.  Input  for  project  name 
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Vessel  Characteristics  obtains  the  basic  dimensions,  form  coefficients, 


LCG,  yaw  radius  of  gyration,  and  bow  profile  information  for  the  hull. 


Figure  30.  Input  for  vessel  characteristics 

Steering  Characteristics  obtains  the  rudder  area,  steering  gear  time 
constant,  position  of  the  rudder  center  of  effort  and  selection  of  single  or  twin  screw. 


Figure  3 1 .  Input  for  steering  characteristics 
Operating  Conditions  obtains  the  water  depth  and  the  initial  vessel  speed. 
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Operating  Conditions 


i 


Figure  32.  Input  for  operating  conditions 

Water  Properties  allows  the  choice  of  fresh  water  @  15C,  salt  water  @ 
15C  or  user-specified  water  properties. 


Figure  33 .  Input  for  water  properties 

2.  Output  from  Maneuvering  Prediction  Program  (MPP1.3) 

When  the  input  process  is  completed,  the  program  can  be  run  through  the 
Analysis  menu  or  by  selecting  the  RunMPP  button.  This  will  produce  a  final  run  input 
window  where  the  selection  is  made  for  the  linear  evaluation  and  there  is  a  need  to  name 
the  run  identifier. 
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Run  MPP  ... 


1 


Options: 

Note: 

W  Linear  Evaluation 

If  you  want  to  check  all  your 

1  Turning  Prediction 

inputs  before  running  MPP. 
please  click  "Verify"  button. 

RudderAngle:  [ 

V  State  Variable  Equations 

degrees 

Verity 

Run  Identifier: 


|  SAW 


Cancel 


OK 


Figure  34.  Input  for  run  identifier 


An  output  report  will  be  generated  as  shown  below.  There  is  a  need  to 
check  that  the  vessel  is  hydrodynamically  open  loop  course  stable  in  order  for  the  values 
of  the  hydrodynamic  derivatives  to  be  usable. 
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Maneuvering  Prediction  Program  -  [i 

MPPl]  I 

[2  File  Edit  View  Input  Analysis  Window 

1  Help 

D|uS  a  J  ®J  #1U 

University  of  Michigan 

Department  of  Naval  Architecture  and  Marine  Engineering 
Maneuvering  Prediction  Program  (MPP-1.3)  by  M.G.  Parsons 


***  Linear  Maneuvering  Criteria  Option  *** 
Reference 


Clarke, D.,  Gedling.P.,  and  Hine,G., 

"The  Application  of  Manoeuvring  Criteria  in  Hull 
Design  using  Linear  Theory,"  Trans.  RINA,  1983 


Run  Identification:  sawl 


Linear  Maneuvering  Derivatives 


Nondimensional  Mass 

M 

prime 

= 

0.017620 

Nondimensional  Mass  Moment 

I 

sub 

zz 

= 

0.001101 

Sway  Velocity  Derivative 

Y 

sub 

V 

= 

-0.025566 

Sway  Acceleration  Derivative 

Y 

sub 

V 

dot 

= 

-0.016110 

Yaw  Velocity  Derivative 

N 

sub 

V 

= 

-0.008402 

Yaw  Acceleration  Derivative 

N 

sub 

V 

dot 

= 

-0.001046 

Sway  Velocity  Derivative 

Y 

sub 

r 

= 

0.005079 

Sway  Acceleration  Derivative 

Y 

sub 

r 

dot 

= 

-0.001200 

Yaw  Velocity  Derivative 

N 

sub 

r 

= 

-0.003691 

Yaw  Acceleration  Derivative 

N 

sub 

r 

dot 

= 

-0.000870 

Sway  Rudder  Derivative 

Y 

sub 

delta 

= 

0.003275 

Yaw  Rudder  Derivative 

N 

sub 

delta 

= 

-0.001605 

Time  Constants  and  Gains  for  Nomoto"s  Equation 


Dominant  Ship  Time  Constant 

T1  prime 

= 

13.7685 

Ship  Time  Constant 

T2  prime 

= 

0.3944 

Numerator  Time  Constant 

T3  prime 

= 

0.8817 

Numerator  Time  Constant 

T4  prime 

= 

0.2791 

1st  Order  Eqn.  Time  Constant 

T  prime 

= 

13.2812 

Rudder  Gain  Factor 

K  prime 

= 

-5.9580 

Rudder  Gain  Factor 

K  sub  v  prime 

= 

3.0506 

Steering  Gear  Time  Constant 

TE  prime 

= 

0.1213 

Evaluation  of  Turning  Ability  and  Stability 


Inverse  Time  Constant 
Inverse  Gain  Factor 


1/|T  prime | 
1/|K  prime] 


0.0753 

0.1678 


Clarke "s  Turning  Index 

Linear  Dynamic  Stability  Criterion 


C  =  0.0000115 
Vessel  is  hydrodynamically  open  loop  course  stable 


Ready 


Figure  35.  Output  from  the  Maneuvering  Prediction  Program  (MPPl  .3) 
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APPENDIX  B.  MATLAB  PROGRAM  CODE 


clear  all 
clc 

o, 

o 

o _ 

O - 

%  Input  Variables 

alpha  =0.3;  %  Increase  in  draft,  values  range  0  to  1 

beta  =0.3;  %  Applied  Fwpb  @  Fwd  area,  values  range  0  to  1 

Delta  =0.;  %  Set  rudder  deflection  angle  to  zero 


%  Declarations 


=  1025.; 

O, 

o 

Density  of  seawater 

9.81; 

o, 

o 

Gravitional  Accelerations 

2  65.  ; 

o, 

o 

Length  of  SAW 

43.2; 

o, 

o 

Breath  of  SAW 

17.38; 

o, 

o 

Draft  of  SAW 

12.86; 

o, 

o 

Velocity  of  SAW  (25knots 

12 . 86m/s) 


O _ 

O - 

%  Formulations 

phi  =  atan  (  (T*alpha) /B) ; 

Fwpb  =  -l*rho*g*L*B*0 . 5* ( (2*T) + (T*alpha) ) *sin (phi) *beta  %  Force 
calculation 

1  =  L/2- ( 0 . 15*L) ;  %  Moment  Length  wrt  amidships 
Nwpb  =  l*Fwpb;  %  Moment  calculation 

Fwl  =  (Fwpb) / (0 . 5*rho* (uA2) * (LA2) ) ;  %non  dimensionize 

Nwl  =  (Nwpb) / (0 . 5*rho* (uA2) * (LA3) ) ;  %non  dimensionize 

Deltal  =  (Delta/180) *pi; 

U  =1.;  %Assume  constant  forward  speed  &  non  dimensionize 

xg  =  0.05; 

o _ 

o - 

%Hydrodynamics  Derivatives  Values  from  Maneuvering  Program  (University 
of  Michigan) 


m 

=  0.017620; 

Iz 

=  0.001101; 

Yv 

=  -0.025566; 

Yvdot 

=  -0.016110; 

Nv 

=  -0.008402; 

Nvdot 

=  -0.001046; 
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Yr  =  0.005079; 

Yrdot  =  -0.001200; 

Nr  =  -0.003691; 

Nrdot  =  -0.000870; 

Ydelta  =  0.003275; 

Ndelta  =  -0.001605; 

Den  =  (Iz-Nrdot) * (m-Yvdot)  -  (m*xg-Yrdot) * (m*xg-Nvdot) ; 

all  =  ( (Iz-Nrdot) *Yv  -  (m*xg-Yrdot) *Nv) /Den; 

al2  =  ( (Iz-Nrdot) * (Yr-m)  -  (m*xg-Yrdot) * (Nr-m*xg) ) /Den; 

a21  =  ( (m-Yvdot) *Nv  -  (m*xg-Nvdot) *Yv) /Den; 

a22  =  ( (m-Yvdot) * (Nr-m*xg)  -  (m*xg-Nvdot) * (Yr-m) ) /Den; 

bl  =  ( (Iz-Nrdot) *Ydelta  -  (m*xg-Yrdot) *Ndelta) /Den; 

b2  =  ( (m-Yvdot) *Ndelta  -  (m*xg-Nvdot) *Ydelta) /Den; 


O _ 

O - 

o, 

o 

%  Initialization 

O, 

o 


v  old 

=  0.; 

r  old 

=  0.; 

psi  old 

=  0.; 

xO  old 

=  0.; 

yO  old 

o. 

=  0.; 

o 

DeltaT 

=  0.001; 

%  Time  step 

increment 

SimTime 

=  5.0; 

%  Simulation 

time 

NT 

=  SimTime/DeltaT; 

%  Number  of 

simulation 

o, 

o 


%  Start  simulation 

o, 

o 

for  i=l :NT, 

O, 

o 

%  External  force 

O, 

o 

Fw  =  Fwl ; 

Nw  =  Nwl ; 

if  (i*DeltaT)  >  2.0 
Fw  =  0; 

Nw  =  0  ; 

end 

if  (i*DeltaT)  <  1.0 
Fw  =  0; 

Nw  =  0; 

end 

cl  =  ( (Iz-Nrdot) *Fw  -  (m*xg-Yrdot) *Nw) /Den; 
c2  =  ( (m-Yvdot) *Nw  -  (m*xg-Nvdot) *Fw) /Den; 

O, 

o 

%  Equations 

O, 

o 

psidot  =  r_old; 
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vdot 

rdot 

xOdot 

yOdot 

o, 

o 

%  First 

O, 

o 

psi_new 

v_new 

r_new 

xO_new 

yO_new 

ang_old 

ang_new 

o, 

o 

%  Store 

o, 

o 


all*v_old  +  al2*r  old  +  bl*Deltal  +  cl; 
a21*v_old  +  a22*r  old  +  b2*Deltal  +  c2; 
(U*cos (psi_old) )  -  (v_old*sin (psi_old) ) ; 
(U*sin (psi_old) )  +(  v_old*cos (psi_old) ) ; 

order  integration  (Euler  -  explicit) 


=  psiold 
=  v_old 
=  r_old 
=  xO  old 
=  yO_old 
=  0; 

=  ang  old 


+  DeltaT*psidot; 

+  DeltaT*vdot; 

+  DeltaT*rdot; 

+  DeltaT*xOdot*u; 

+  DeltaT*yOdot*u; 

(atan (abs (y0_new) /xO^new) ) *57.2958; 


results 


psi (i) 

=  psi  new; 

psidot  v(i) 

=  psidot; 

v  (i) 

=  v  new; 

vdot  v(i) 

=  vdot; 

r  (i) 

=  r  new; 

rdot  v(i) 

=  rdot; 

xO  (i) 

=  xO  new; 

y0  (i) 

=  yO  new; 

Fw  v ( i ) 

=  Fw; 

Nw  v ( i ) 

=  Nw; 

time ( i ) 

=  i*DeltaT; 

angle  ( i ) 

=  ang  new; 

o, 

o 


%  Update  state  vector  x  for  the  next  simulation  step 

o, 

o 


psi_old 

v_old 

r_old 

x0_old 

y0_old 

end 


psi  new; 
v  new; 
r  new; 
xO  new; 
y0_new; 


O, 

o 

%  Plotting  of  graphs 

O, 

o 


figure ( 1 ) 

plot (time, psi*57.2958) ; xlabel ( ' t  ( sec) ' ) ; ylabel ( ' \psi  (deg) ' ) ; grid 

o, 

o 


figure  (2 ) 

plot (time, v) ; xlabel ( ' t  ( sec)  ' ) ; ylabel ( ' v ' ) ; grid 

o, 

o 


figure ( 3 ) 

plot (time, r ) ; xlabel ( ' t  ( sec)  ' ) ; ylabel ( ' r ' ) ; grid 

O, 

o 

figure ( 4 ) 

plot (time, psidot_v) ; xlabel ( ' t  ( sec)  ' ) ; ylabel ( '  (d/dt ) \psi ' ) ; grid 

o, 

o 

figure ( 5 ) 

plot (time, vdot_v) ; xlabel ( ' t  ( sec)  ' ) ; ylabel ( '  (d/dt ) v '); grid 
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figure ( 6 ) 

plot (time, rdot_v) ; xlabel ( ' t  (sec) ' ) ; ylabel ( ' (d/dt ) r ' ) ; grid 

o, 

o 

figure ( 7 ) 

plot (time, Fw_v) ; xlabel ( ' t  (sec) ' ) ; ylabel ( ' Fw ' ) ; grid 

O, 

o 

figure ( 8 ) 

plot (time, Nw_v) ; xlabel ( ' t  ( sec) ' ) ; ylabel ( ' Nw ' ) ; grid 

o, 

o 

figure ( 9 ) 

plot (xO , yO ) ; xlabel ( ' x ' ) ; ylabel ( ' y ' ) ; grid 
title  ('xO  vs  yO  (Dimensional) ') 

%legend (' Force  #1', 'Force  #2', 'Force  #3', 'Force  #4 ',' location ',' best ' ) 

o, 

o 

figure ( 10 ) 

plot (time, angle,  ' b- ' ) ; xlabel ( ' t  ( sec)  ' ) ; ylabel ( ' deg ' ) ; grid 
title  ('Time  vs  Angle  of  Deflection  (Dimensional) ') 

%legend (' Force  #1',  'Force  #2',  'Force  #3',  'Force  #4 ',' location ',' best ' ) 
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